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Os| ■ Abstract 

Using a sample of 9.7 x 10 6 BB meson pairs collected with the CLEO detector, 

p^j \ we study inclusive B meson decays to the Xd and Xc2 charmonia states. We 

measure the branching fraction for the inclusive Xd production in B decays to 
u* n(Tz _>. •,, , v\ — (a i/i -u n qi -t- n a(\\ ^ m-3 



> 



o 



X 



be B(B -> XdX) = (4.14 ± 0.31 ± 0.40) x 10 -3 , where the first uncertainty is 



^D ■ statistical and the second one is systematic. We obtain the branching fractions 



for direct Xc\ and Xc2 production in B decays by subtracting the contribution 
^ \ from the decay chain B — > ip(2S)X with *p(2S) — » Xci,27- We measure B{B — » 

Xci[direct]X) = (3.83±0.31±0.40)xl0~ 3 . No statistically significant signal for 
Xc2 production is observed in either case. We determine the 95% C.L. upper 
limits to be B(B -» Xc2X) < 2.0 x 10" 3 , B(B -» Xc2[direct]X) < 1.7 x 10~ 3 , 
D ! and B(B -> ^ [direct] X)/£(.B -» Xci[direct]X) < 0.44. All quoted results 

Oh 



Qh' are preliminary. 
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The recent measurements of charmonium production in various high-energy physics reac- 
tions have brought welcome surprises and challenged our understanding both of heavy-quark 
production and of quarkonium bound state formation. The measurement of a large produc- 
tion rate of high-Pr charmonium at the Tevatron [|l[] was in sharp disagreement with the 
then-standard color-singlet model. The development of the NRQCD factorization frame- 
work [[J has put the calculations of the inclusive charmonium production on a rigorous 
footing. The high-P^ charmonium production rate at the Tevatron is now well understood 
in this formalism. The recent CDF measurement of charmonium polarization ||, however, 
appears to disagree with the NRQCD prediction. The older color-evaporation model accom- 
modates both the high-Pr charmonium production rate and polarization measurements at 
the Tevatron Q. 

Inclusive B meson decays to charmonia is another area to confront theoretical predic- 
tions with experimental data. The color-singlet contribution, for example, is a factor of 5-10 
below the observed inclusive J /if) production rate . A measurement of the Xc2-to-Xci 
production ratio in B decays provides an especially clean test of the charmonium produc- 
tion models. The V — A current C7 M (1 — 75)0 cannot create a cc pair in a 2S+1 Lj = 3 P 2 
state, therefore the decay B —> Xc2X is forbidden at leading order in a s in the color-singlet 
model @ . The importance of the color-octet mechanism for the \ c production in B decays 
was recognized || even before the development of the NRQCD framework. The NRQCD 
calculations of the B decays to charmonia at the next-to-leading order in a s can be found 
in Ref. |J. The color-octet contribution to B — > XcjX decays is proportional to 2.1 + 1. 
Therefore, if this mechanism dominates, then the Xc2-to-Xd production ratio should be 5:3. 
On the contrary, if the color-singlet contribution dominates, then the ratio should be 0:1. 
The color-evaporation model predicts Xc2 '■ Xci = 5:3 ||. 

Our data were collected at the Cornell Electron Storage Ring (CESR) with two configu- 
rations of the CLEO detector called CLEO II and CLEO II. V 0. The components of 
the CLEO detector most relevant to this analysis are the charged particle tracking system, 
the Csl electromagnetic calorimeter, and the muon chambers. In CLEO II the momenta 
of charged particles are measured in a tracking system consisting of a 6-layer straw tube 
chamber, a 10-layer precision drift chamber, and a 51-layer main drift chamber, all operat- 
ing inside a 1.5 T solenoidal magnet. The main drift chamber also provides a measurement 
of the specific ionization, dE/dx, used for particle identification. For CLEO II. V, the straw 
tube chamber was replaced with a 3-layer silicon vertex detector, and the gas in the main 
drift chamber was changed from an argon-ethane to a helium-propane mixture. The muon 
chambers consist of proportional counters placed at increasing depth in the steel absorber. 

We use 9.2 fb" 1 of e + e" data taken at the T(45) resonance and 4.6 ftr 1 taken 60 MeV 
below the T(4S) resonance (off-T(45 l ) sample). Two thirds of the data were collected with 
the CLEO II. V detector. The simulated event samples used in this analysis were generated 
with a GEANT-based [|12] simulation of the CLEO detector response and were processed in 
a similar manner as the data. 

We reconstruct the xa and Xc2 radiative decays to J/ if). The branching fractions for the 



Xci,2 —> J/ipl decays are, respectively, (27.3 ± 1.6)% and (13.5 ± 1.1)% [jD|, whereas the 
branching fraction for the Xco —> J/ipl decay is only (0.66 ± 0.18)%. In addition, the Xco 
production rate in B decays is expected to be smaller than the x&,2 rates [|]||. We therefore 
do not attempt to measure Xco production in this analysis. 



We reconstruct both J/tp — > e + e~ and J/ip — ► yU + /i~ decays. Electron candidates are 
identified based on the ratio of the track momentum to the associated shower energy in 
the Csl calorimeter and on the dE/dx measurement. The internal bremsstrahlung in the 
J/ip — > e + e~ decay as well as the bremsstrahlung in the detector material produces a long 
radiative tail in the e + e~ invariant mass distribution and impedes efficient J/ip — > e + e~ 
detection. We recover some of the bremsstrahlung photons by selecting the photon shower 
with the smallest opening angle with respect to the direction of the e ± track evaluated at the 
interaction point, and then requiring this opening angle to be smaller than 5°. We therefore 
refer to the e + ('j)e~('~f) invariant mass when we describe the J/ip — ► e + e~ reconstruction. 
For the J/ip — > /U + /i~ reconstruction, one of the muon candidates is required to penetrate 
the steel absorber to a depth greater than 3 nuclear interaction lengths. We relax the 
absorber penetration requirement for the second muon candidate if it is not expected to 
reach a muon chamber either because its energy is too low or because it does not point 
to a region of the detector covered by the muon chambers. For these muon candidates we 
require the ionization signature in the Csl calorimeter to be consistent with that of a muon. 
We use the normalized invariant mass for the J/ip candidate selection. For example, the 
normalized J/ip — > jjl + h~ mass is defined as [M(/i + /z~) — Mj/^]/a(M), where Mj/^ is the 



world average value of the J/ip mass |T3| and a(M) is the calculated mass resolution for that 
particular n + fi~ combination. The average £ + £~ invariant mass resolution is 12 MeV/c 2 . 
The normalized mass distributions for the J/ip — > £ + £~ candidates are shown in Fig. |l|. We 
require the normalized mass to be from —6 to +3 for the J/ip — > e + e~ and from —4 to +3 
for the J/ip — > /U + /i~ candidates. The momentum of the J/ip candidates is required to be 
less than 2 GeV/c, which is slightly above the maximal J/ip momentum in B decays. 
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FIG. 1. Normalized invariant mass of the (a) J/ip — > e + e~ and (b) J/ip — > /ur [i~ candidates in 
the data. The momentum of the J/ip candidates is required to be less than 2 GeV/c. The shaded 
histogram represents the luminosity-scaled data taken 60 MeV below the T(4S) showing the level 
of background from non-BB events. 



Photon candidates for Xa,2 — > J/ipl reconstruction must be detected in the barrel region 
of the calorimeter, defined as | cos 6^1 < 0.71, where 6* 7 is the angle between the beam axis 



and the candidate photon. Most of the photons in T(45 r ) — ► BB events come from tt° decays. 
We therefore do not use a photon for the Xa,2 —* J/ipl reconstruction if it can be paired 
with another photon to produce a tt° candidate with the normalized 7r° — * 77 mass between 
-3 and +2. 

We determine the Xci and Xc2 yields in a maximum-likelihood fit to the mass- difference 
distribution M{J/if) r y) — M(J/ip) (Fig. ||), where M(J/ip) is the measured mass of a J /if) 
candidate. The excellent electromagnetic calorimeter allows us to resolve the xa and Xc2 
peaks. The M ( J/i/jj) — M(J/ip) mass-difference resolution is approximately 8 MeV/c 2 and 
is dominated by the photon energy resolution. The background in the fit is approximated 
with a 5th order Chebyshev polynomial; all the polynomial coefficients are allowed to float. 
The Xci and Xc2 signal shapes are fit with templates extracted from Monte Carlo simulation; 
only the template normalizations are floating in the fit. The Xci and Xc2 signal yields in 
the T(4S) data are iV ON (x c i) = 672 ± 47[stat] and iV ON (x c2 ) = 83 ± 37[stat]. The xa 
and Xc2 yields in off-T(45) data are consistent with zero: iV OFF (x c i) = 4.1 ± 7.1[stat] and 
iV OFF (x C 2) = 1.1 ±6.5[stat]. Subtracting the contributions from noia-BB continuum events, 
we obtain the total inclusive B — > XaX and B — > Xc2X yields 

N(B -► XdX) = 664 ± 49[stat] and N(B -» Xc2*) = 81 ± 39[stat]. 
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FIG. 2. Plot (a) shows the M(J/i/)j) - M(J/i/j) distribution in the T(45) data (points with 
error bars). The fit function is shown by a solid line with the background component represented 
by a dashed line. Plot (b) shows the background-subtracted distribution in the signal region with 
the Xci and Xc2 fit components represented by a solid line. 



Taking into account the systematic uncertainties associated with the fit, we determine the 
B — > Xc2X signal yield significance to be 2.0 standard deviations. Subtracting the feeddown 
from the decay chain B — ► if>(2S)X with ip(2S) — > Xc2l and accounting for the associated 
systematic uncertainty, we likewise determine the significance of the evidence for the direct 
B —* Xc2X decays to be 1.4 standard deviations. 



To calculate branching fractions for the B — » x&pX decays, we use the measured signal 
yields N(B — > Xci^X), reconstruction efficiencies extracted from simulation, the number of 
produced BB pairs, and the daughter branching fractions. For the calculation of the rates 
for direct B — > Xci,2X decays, we make an assumption that the only other source of Xci,2 
meson production in B decays is the decay chain B — > ip{2S)X with ip{2S) — > Xd,2l- The 
resulting branching fractions are listed in Table |. 

TABLE I. The measured branching fractions for inclusive B decays to Xci an d Xc2- We subtract 
the ip(2S) — > X c i27 feeddown to arrive at the branching fractions for direct B — > Xci,2X decays. 
The correlation between the uncertainties is taken into account in the calculation of the Xc2-to~Xci 
production ratio (last row). 



Branching fraction 
or ratio 


Measured value 


95% C.L. 
upper limit [15 


B(B -» XdX) 
B(B ^ X ci[direct}X) 


(4.14 ±0.31 ±0.40) x 10~ 3 
(3.83 ±0.31 ±0.40) x 10~ 3 


Z 


B(B -» Xc2X) 
B(B ^ X c2[dhect]X) 


(0.98 ±0.48 ±0.15) x 10" 3 
(0.71 ±0.48 ±0.16) x 10~ 3 


< 2.0 x 10" 3 a 

< 1.7 x 10~ 3 


B(B ^ X c2[dircct\X) 
B(B ^ X ci[dircct]X) 


0.18 ±0.13 ±0.04 


<0.44 



a The 95% confidence interval in the "unified" approach G5J is [0.2; 2.0] x 10 
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The systematic uncertainties are listed in Table |T|. The sources of the uncertainty can 
be grouped into three categories: 

(i) Fit procedure. — This category includes the systematic uncertainties due to the choice 
of the signal and background shapes as well as the bin size. To fit the Xci and Xc2 signal, we 
use the templates extracted from simulation. We therefore are sensitive to the imperfections 
in the simulation of the photon energy measurement. The systematic uncertainties associated 
with the simulation of the calorimeter response are estimated by comparing the 7r° — > 77 
invariant mass lineshapes for inclusive 7r° candidates in data and Monte Carlo samples. Then 
we modify the x& an d Xc2 templates to determine the resulting uncertainty in the signal 
yields. To estimate the uncertainty associated with the calorimeter energy scale, we shift the 
Xci and Xc2 templates by ±0.6 MeV/c 2 in the fit. The uncertainty due to time-dependent 
variations of the calorimeter energy scale is small compared to the overall energy scale 
uncertainty. To estimate the uncertainty due the calorimeter energy resolution, we change the 
width of the Xci an d Xc2 templates by ±4%. The uncertainty in background shape is probed 
by fitting the background with the template extracted from the high-statistics samples of 
simulated T(45 r ) — * BB and non-BB continuum events; only the template normalization, 
not its shape, is allowed to float in the fit. 

(ii) Efficiency calculation. — This category includes the uncertainties in the number of 
produced BB pairs, tracking efficiency, photon detection efficiency, lepton detection effi- 
ciency, statistical uncertainties of the simulated event samples, and the model-dependence 
in the simulation of B — > Xci^X decays. The angular distribution of the Xc\,2 -^ J/^l decays 
affects the photon energy spectrum. We define the helicity angle 6 h to be the angle between 
the 7 direction in Xc frame and the Xc direction in the B frame. We assume flat cos 9h dis- 

7 



tribution in our simulation of B — > Xci,2^ decays. Systematic uncertainty associated with 
this assumption is estimated by comparing the reconstruction efficiencies in the Monte Carlo 
samples with I(9h) oc sin 2 6h and I(9h) oc cos 2 9h angular distributions. Another source of 
uncertainty is the modeling of the X system in the simulation of B — > Xa^X decays. The 
photon detection efficiency depends on the assumed model through the Xc momentum spec- 
trum and the ir° multiplicity of the final state. In our simulation, we assume that X is 
either K or one of the higher K resonances: 24% K, 24% Hf (892), 14% 7^(1270), 14% 
#i(1400), 13% K*(1430), and 11% K*(U30); we also include the decay chain B -> *p(2S)X 
with '0(25') — > Xci,2l- To estimate the systematic uncertainty, we compare the Xc —> J /''Pi 
detection efficiency extracted using this sample with the efficiency in the sample where we 
assume that X is either a K + or Kg — > 7r + 7r~. 

(iii) Assumed branching fractions. — - This category includes the uncertainties on the 
branching fractions. We use the following values of the daughter branching fractions: 
B(J/i/j -> £+£-) = (5.894 ± 0.086)% fTJ, B( X d -> J/i>i) = (27.3 ± 1.6)% fTJ, and 
B{Xc2 -► J/iH) = (13.5 ± 1.1)% 0. In calculation of B[B -► X d,2 [direct] X) we also 
assume the following values: B(S -> ij}{2S)X) = (3.5 ±0.5) x 10~ 3 [[TJ, i3(^(25) -+ Xal) = 
(8.7 ± 0.8)% H, and 6(^(25) -► Xc27) = (7.8 ± 0.8)% |3|. 



TABLE II. Systematic uncertainties on the measured branching fractions. 



Source of 


relative uncertainty in 


% 


systematic uncertainty 


£(B - Xci*) 


B(B 


- Xc2X) 


(i) Fit procedure 








Calorimeter energy scale 


0.4 




5.6 


Calorimeter energy resolution 


2.8 




6.9 


Background shape 


1.8 




6.8 


Bin size 


0.0 




1.9 


(ii) Efficiency calculation 








N(BB) 


2.0 




2.0 


Tracking efficiency 


2.0 




2.0 


Lepton identification 


4.2 




4.2 


Photon finding 


2.5 




2.5 


Monte Carlo statistics 


0.7 




0.7 


Composition of X in B — > Xci^X simulation 


3.3 




3.3 


Angular distribution for Xci,2 — * J /in 


1.0 




1.0 


(iii) Assumed branching fractions 








B(Xci,2 - W?) 


5.9 




8.1 


B(J/V> -» ^~) 


1.5 




1.5 


B(B -► V(25)A) a 


1.1 




5.5 


£(V(2S) -> Xc i,27) a 


0.7 




4.0 



a Contributes only to uncertainty on B{B — ► Xci,2 [direct] X). 



In conclusion, we have measured the branching fractions for inclusive B meson decays to 
Xci and Xc2 charmonia states. Our measurements are consistent with and improve upon the 
previous CLEO results B. Our measurement of the branching ratio for direct Xc2 and xa 



production in B decays is consistent with the prediction of the color-singlet model J/] and 
disagrees with the color-evaporation model ||. In NRQCD framework, our measurement 
suggests that the color-octet mechanism does not dominate in B — > \cX decays. 
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